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ABSTRACT 

The spectral properties of the Dirac Hamiltonian in the the Aharonov-Bohm potential are 
discussed. By using the Krein-Friedel formula, the density of states (DOS) for different 
self-adjoint extensions is calculated. As in the nonrelativistic case, whenever a bound state 
is present in the spectrum it is always accompanied by a (anti)resonance at the energy. 
The Aharonov-Casher theorem must be corrected for singular field configurations. There 
are no zero (threshold) modes in the Aharonov-Bohm potential. For our choice of the 2d 
Dirac Hamiltonian, the phase-shift flip is shown to occur at only positive energies. This 
flip gives rise to a surplus of the DOS at the lower threshold coming entirely from the 
continuous part of the spectrum. The results are applied to several physical quantities: 
the total energy, induced fermion-number, and the axial anomaly. 
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1. Introduction.- In this paper the change Ap a (E) over all space of the density of 
states (DOS) induced by the Aharonov-Bohm (AB) potential A(r) [1], 

A = 0, A v = ^- = ^§ , (1) 
Zirr Zirr 

in the radial gauge, is calculated for the two-dimensional (2d) Dirac Hamiltonian. The 
AB potential can be uderstood in a more general sense here. Usually $ = a $ is the 
total flux through the flux tube and $ is the flux quantum, $ = hc/\e\. However, the 
same potential (1) is induced by a cosmic string provided the identification a = e/Qniggs 
is made, where e and Qiiiggs are respectively the charge of a test particle and the charge 
of the Higgs particle [2, 3]. In what follows, n and r\ will be respectively the integer and 
the fractional part of a = n + r\. 

If one neglects anomalous magnetic moment of the electron, the Dirac Hamiltonian in 
the AB potential (and in any magnetic field) considered as an operator in L 2 is a formal 
operator with supersymmetry [4]. A contradiction between periodicity of the spectrum 
under the substitution a — > ail (as long as a does not change its sign) and the Aharonov- 
Casher (AC) theorem [5] is resolved. One finds that there are no threshold (zero) modes 
in the AB potential for any a (see Fig. 1). This can be confirmed independently by 
the calculation of the Witten index [6]. Similarly as in the nonrelativistic case [7], our 
calculation of the contribution of continuous spectrum to the integrated DOS (IDOS) 
will be based on the analysis of phase shifts and use of the Krein-Friedel formula [8, 9] 
regularized by the ('-function [7]. First, we shall consider a generic self-adjoint extension, 
when the supersymmetry of the Dirac Hamiltonian is lost and a bound state with energy 
Eb occurs in the spectrum in the channel I = —n — 1. We shall calculate the change of the 
IDOS and show that the bound state is always accompanied by a (anti)resonance [10]. In 
the supersymmetric case, after taking the limits Eb —> ±m in the channel I = —n — 1, the 
bound state disappear and instead of the resonance the phase-shift flip (cf. [11]) occurs. 

The DOS provides an important link between different physical quantities. The knowl- 
edge of Ap a (E) determines the spectral asymmetry a a (S), 

a a {£) = Ap a {£) - A Pa {-£), (2) 

where £ = \E\ > 0. In contrast to three dimensions, the spectrum of the 2d Dirac Hamil- 
tonian is in general asymmetric and, if supersymmetry holds, the spectral asymmetry is 
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Figure 1: In the case of a 2d Dirac Hamiltonian in a generic finite-flux magnetic field the 
threshold states occur at one and only one threshold. There are, however, no threshold 
(zero) modes in the AB potential. In the latter case, the spectral asymmetry is entirely 
accounted for by the continuous part of the spectrum. The phase-shift flip which occurs 
at positive energies in the I = —n — 1 channel gives rise to the net surplus of 77 states at 
the lower threshold. 



concentrated at the thresholds. The spectral asymmetry measures how the vacuum is 
distorted under the presence of an external field, and therefore it is natural that spectral 
asymmetry cr a (£) determines the one-loop contribution to the effective energy, the 
induced fermion number Q, and, in the case of the massless Euclidean Dirac operator, the 
axial anomaly A. One has 



poo poo poo 

E 1 eff = J o £a a {S)dS, Q = -^J o v a {£)dS, A = j o a a {S)dS, 



(3) 



and Q = —A/2 [12, 13]. By using the above relations one can check for the consistency 
of the result for either of these quantities. 

2. The Dirac Hamiltonian in the AB potential- The Dirac Hamiltonian in a magnetic 
field (in units with h~ = c = 1) is 



H(A) 



D —m 



D = E(?i + A + 2 (?2 + ^2), 



(4) 



where pj = idj. We have taken the standard representation of 7 matrices in 2d, 7 = (3 = 
cr 3 , 7 1 = ia 2 , and = —icri, where aj are the Pauli matrices. In what follows, the charge 
e will be assumed to be negative, e = — |e|. 

In the AB potential, after separation of variables in polar coordinates, H(A) reduces 
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to the direct sum, H(A) = ©j/i m ,j, of channel operators h m j in L 2 [(0, oo), rdr], 



m —i [d r + 

-i ld r — - ) —m 



(5) 



where v = I + a [3]. One finds immediately that h* m ^m^-m = —h mi i, where '*' means the 
complex conjugation. Therefore, for E = —£ < — m, the scattering states are given by 

r >¥ >) = r >¥ >)| m _>_ m . (6) 

2.1. Point spectrum and the Aharonov-Casher theorem.- Formally, the set (5) of radial 
Hamiltonians is invariant under the substitution a — > a± 1. As long as a does not change 
its sign, the spectrum of H(A) must be periodic in a (see below). This is, however, 
seemingly in contradiction with the AC theorem [5] which states that the point spectrum 
is a function of a. In a general finite-flux magnetic field B(r), 

[ B(r)d 2 r = $ = const, (7) 

the AC theorem tells us that H(A) has exactly either a — 1 = n — lorn threshold states at 
one and only one of the thresholds E = ±m (in the present case at E = —m) depending 
whether the flux a is an integer or not. In the case of the 2d massless Euclidean Dirac 
operator, the threshold states are actually zero modes. 

To resolve this contradiction, one must check for validity of the AC theorem in the 
case of a singular flux tube with zero radius. In the latter case, one must check for square 
integrability not only at infinity but at the position of singularities of the field, too. It 
is here where the theorems fail. If one takes proof of the AC theorem as presented, for 
example in Ref. [4], p. 198, and substitutes for B(r), 



B{r) = 2Tra5( 2 \r) = -5{r), (8) 

r 

one finds immediately that threshold (zero) modes are 



*(') *-*W( Xl - ix2 ) } . . . , e -*W( Xl - ix 2 ) n ~\ (9) 



e x ' , e 



where 



0(r) = ^- j^\n\r -r'\B{r')d 2 r' = aln\r\. (10) 
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Threshold (zero) modes are obviously singular at the origin and are not elements of 
L 2 [(0, oo), rdr]. One can show this directly by using h mt i as well. For a threshold state at 
E = m to exist, the upper component xi °f the Dirac spinor has to obey 

[d r - xi(r) = 0. (11) 

At the threshold E = —m one then obtains 

(d r + ^) Xa (r) = 0. (12) 

These two equations can be easily integrated. The solutions are 

Xi{r)=r" and Xa (r) = , (13) 

and they are obviously not in L 2 [(0, oo), rdr] for any I. If they are square integrable at 
the infinity they are not so at the origin and vice versa. 

2.2. Continuous spectrum.- In order to use the Krein-Friedel formula [8] for the 
continuous part of the spectrum one must determine phase shifts. Scattering states 
^£,1 = x{ r ) e%l ' f ' e ~ %£t ^ h f° r E > m are given in terms of Bessel functions with 

X[r) N { zVS^l (e,)^ 1 J ei(l/+1) (fcr)e* J ' [ " V 

N being a normalization factor, and £/ = ±1 [3]. The sign of £/, except for the channel 
I = —n — 1, is fixed by the square integrability condition at the origin [3]. The two- 
component solutions of the massive Dirac equation have only one degree of freedom that 
is reflected in the equality of up and down phase shifts [14]. Except for the channel 
I = —n — 1, phase shifts are invariant under the change of the sign of the energy, E — > —E, 
and are given by 

5? = 5?=\ -™> \ > ~ n ~ l : (15) 
' ' [ ira, I < -n - 1. K 1 

In the channel I = —n — 1 two different solutions are possible, 

X-(r) = -( ^ £ + m J n-^r) ) (lg) 



anc 



X + (r) = -( >/ f + m J i->i( kr ) . ) (17) 



with opposite phase shifts, 

S~ = ~^V, <$ + = ^V- ( 18 ) 

Note that state % + (x~) 1S regular (singular) in the limit 77 — >• + and singular (regular) in 
the limit 77 — > 1_. Provided that one takes solution % + for I = —n — 1, then phase shifts 
for all I are given by Si = ^ir(\l \ — \l -\- a\). This ambiguity in the channel I = —n — 1 is the 
consequence of the fact that the AB potential is in the so-called limit circle case at the 
origin (see Ref. [15], p. 152) where the boundary conditions have to be specified. A given 
combination of % + and %~ then characterizes a particular self-adjoint extension of H(A) 
[3]. Provided that a = —n — r\ < 0, the critical channel is I = n. Therefore, despite the 
fact that the set (5) of radial Hamiltonians is invariant under the substitution a —> ail, 
the spectrum is not. The spectrum is periodic only as long as a does not change its sign. 

3. Self-adjoint extensions.- Distinguished feature of a generic self-adjoint extension is 
that the point spectrum is not empty and a bound state can occur with energy inside the 
mass gap (— m,m) [3, 16]. Both supersymmetry and the scale invariance are broken in 
this case. The bound state of energy Ef, in the channel I = —n — 1 is given in terms of 
modified Bessel functions K u (z) [17], 

where k = \Jm 2 — E£. In the presence of a bound state of energy E{,, scattering states 
with energy E > m take a general form 

V S; - n -i{t,r,<p) = [ X + {r) sin fi + X -{r) cos fi] e -^+iW e ^st/n^ (2Q) 

where x + ( r ) an d X~( r ) are given respectively by (16) and (17). In contrast to [3], we 
shall not use the factor ( — l) n+1 for x~( r )- Such as in the non-relativistic case [7], in the 
presence of a bound state of energy Ef, the parameter \i of the self-adjoint extension is a 
function of energy E, 

tan u = — h-^-z -• 21 

P {m 2 -ElY £-m K 1 

The formula is only valid for E > m. Scattering states at energies E < —m are obtained 

directly from positive energy states by using relation (6). If necessary, fi± will be written 



to distinguish between positive and negative energies. A close inspection of (21) reveals 
that under m — > —m both the sign and the value of tan fi + change, 

m + E b (S 2 - m 2 f . . 

tanu+ — >■ tanu_ = —- — — rr (22) 

P+ P (m 2 - ElY £ + m y 1 

The change amounts to the following two changes in (21), E b —> —Eb, and £ — > —£. The 
corresponding phase shift is given by 

S-n-i(E) = l -ita + A.^E), (23) 

where 

sin(r77r) 

A_ n _i = - arctan (24) 

cos(r77r) — tan \i 

The analysis presented here again shows consistency of our result about the absence 
of threshold (zero) modes in the AB potential. As we have seen, the I = —n — 1 channel is 
the only channel where the sign of factor £/ in (14) is not fixed by the square integrability 
condition at the origin and where, in principle, a bound state can appear. On the other 
hand, the AC theorem allows for zero modes only in the channels I = 0, —1, ... , —n + 1 
[cf. (9)], none of which is the critical channel I = —n — 1. 

4. The density of states and the Krein-Friedel formula.- Let us denote by Ny(E) and 
No(E) respectively the IDOS induced by the presence of a scatterer V and the free IDOS. 
The contribution of scattering states to the change ANv(E) of the IDOS induced by the 
presence of a scatterer with a finite range of interaction is known to be determined by the 
Krein-Friedel formula [8] directly as the sum over all phase shifts, 

AN V (E) = N V {E) - N (E) = - V 5 t {E) = {2-ki)- 1 In det S, (25) 

* i 

S being the total on-shell S-matrix. In a recent papers [7] we have shown that the Krein- 
Friedel formula can have a wider range of applicability. In particular, we have established 
that the Krein-Friedel formula can also be used for the long-ranged AB potential, provided 
that the sum over phase shifts is regularized by the ('-function. The fact that phase shifts 
can be rather easily calculated without any care of the proper normalization of wave 
functions greatly facilitates the calculation. Moreover, by means of the Krein-Friedel 
formula it is rather easy to calculate the change of the IDOS for all possible self-adjoint 
extensions of (4). In the AB potential, the IDOS N a (E) is as usually defined by N a (E) = 



/foo p a (E') dE' , where p a (E) = — (l/7r)ImTrG ! a (x,x, E + ie), and G a (x,x, £ + ie) is the 
resolvent of H(A) in the AB potential with the flux corresponding to a. By using the 
('-function regularization [7] one finds in each energy continuum that 

oo oo 

lndetS= J2 2iS i = 2iA -n-i + ™ J2 ~ \ l + a \) 

l= — oo l= — oo 

oo oo oo 

= 2lA_ n _ 1 + 27T 2 £ Z" S - £(Z + 7?)- S - £(Z - 7/)- 

. z=i z=o z=i 

= 22A_ n _! + 27T [2Cfl(s) " Ch(s, V ) ~ Ch(s, 1 - v)]\.=-l 

= 2iA_ n _ 1 -nr V (l-r 1 ), (26) 

where (j? and (h are the Riemann and the Hurwitz (-functions. Therefore, according to 
the Krein-Friedel formula (25), 

sin(777r) 



AN a (E) = -\-q{\ -77) - -arctan- 

z 7T cos (7771 J — tan fj, 



One can repeat the calculation for a < to verify that 

AM w (-£) = AW W (£). 



(27) 



(28) 



4.1. TTie resonance.- tan is positive for E > m [see (21)] and hence there is a 
typical (anti)resonance in the relativistic case for < 77 < 1/2 at the energy £ res which 
satisfies 

f -m 1 m-E h (2Q) 



m 



m 2 Y 



COS 777T (m* 



£ fc 2 )" 



The phase shift 5- n -i(E) (23) changes by — it in the direction of increasing energy and 
the integrated DOS (27) has a sharp decrease by one. For 1/2 < 77 < 1 the cosine in 
(27) changes its sign and the antiresonance disappears at positive energies [18]. Instead, 
a resonance occurs at the negative energy continuum at the energy E res = —£ res which is 
determined by the equation 

£res + rn _ 1 m + E b 
(£ 2 es - m 2 )" " I cos 777r| (m 2 - £ 2 )" ' [ ' 

77 = 1/2 is a special point since the (anti)resonance is at infinity. Obviously, the shape of 
the (anti)resonance is not of the Breit-Wigner type. 

4.2. Supersymmetric case and the phase- shift flip.- One knows that in the supersym- 
metric case there is neither bound state nor the (anti)resonance. Therefore, to get rid of 
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the bound state we shall consider limits Eb —> ±m. If Eb changes, so do E res and tan//±. 
Limiting values of E res and tan p,± are presented in Tab. I. 



TABLE I. Limiting values of E res and tan p,± when Eb —> ±m. 





< 77 < 1/2 


1/2 < 77 < 1 


< < 1/2 


1/2 < 77 < 1 


E b t m 


E res -j. TTL 




tan //_)_—)■ 


tan — >■ —00 


Eb I —m 


E res | oo 




tan fJ-+ —> oo 


tan (J-- —> 



Formally, both cases £^ ^ m and £k ^ —m seem to be plausible. In each case, scatter- 
ing states in one of the energy continuums are described by state % + and in the other 
continuum by state %~. As a result, one has the phase-shift flip in one energy continuum, 

5 + ->■ 5~ = ~^V, 

and the phase shift 5 + in the second energy continuum. In order that our results be 
compatible with previous calculations of the spectral aymmetry and the axial anomaly 
[13], the limit Eb t m was chosen. In this case, the phase-shift flip occurs at positive 
energy continuum while at the negative energy continuum the phase shift is 5 + = irr]/2. 
Therefore, the contribution of scattering states with energy E > m to the change of the 
DOS is 

A Pa (E) = Pa {E) - p (E) = -^77(3 - 77) 5{E - m). (31) 
while the contribution of scattering states with E < —m to the change of the DOS is 

Ap a (E) = -^77(1 - V )5(E + m). (32) 

As expected, the change of the DOS is concentrated at the thresholds where it is propor- 
tional to delta functions. 

5. Spectral asymmetry , effective energy, induced fermion number, and the axial anoma- 
ly.- In the absence of the bound state, when supersymmetry holds, the spectral asymmetry 
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is entirely determined by the continuous part of the spectrum. One finds from (31) and 
(32) that the contribution is actually determined by the phase-shift flip, 

a a (£) = -r]S(£-m). (33) 

Correspondingly, according to (3), 

E^^-mrj, Q = \v, A=-rj. (34) 

The phase-shift flip which occurs at positive energies in the I = —n — 1 channel gives 
rise to the net surplus of r\ states at the lower threshold with respect to the upper energy 
threshold (see Fig. 1). When a reaches 1, then, in the limit a — > 1_, there is seemingly a 
surplus of one state at the lower threshold with respect to the upper threshold. However, 
for a = 1 all quantities of interest, including basic set (5) of radial Hamiltonians and 
wave functions (14) are identical to those in the case for a = 0. By using the fact that 
phase shifts are only defined modulo 7r, one finds that phase shifts depend periodically 
(although discontinuously) on a [cf. (15)]. Therefore, we argue that spectral asymmetry 
a a , and consequently E^, Q, and A change discontinuously to their original zero values 
for a = 1. This can already be seen from our results (31-32) for the change of the DOS. 
Actually, we have already encountered a discontinuous behaviour in the discussion of the 
phase-shift flip. Despite that %~ is singular in the limit 77 — >- 0, it is this state which occurs 
in the spectrum for E > m, and not state % + . 

Another justification for our conclusion for the discontinuous behaviour of the spectral 
asymmetry comes from the nonrelativistic case. In the latter case one knows that there 
exists one physical quantity closely associated with the DOS, and hence with the spectral 
asymmetry, which changes discontinuously with the flux a [19]. This experimentally 
measurable quantity [20] is the so-called persistent current induced by the AB potential 
[21]. 

In the case for a generic self-adjoint extension when the bound state (19) is present, 
the supersymmetry is lost, because spectral asymmetry a a is no longer concentrated at 
the mass thresholds, and respective restrictions of D on (Ker D~) L and on (Ker D^) 1 - 
are no longer unitary equivalent. By using results (27) for the IDOS one finds 

a a (S) = sign(£: fc ) 5(8 - \E b \) + ~ { [AV^) - AI^^)] e(S - m)} , (35) 

9 



where is the Heaviside step function, and A_ n _ 1 (S) is the change of the conventional 
phase shift (24) respectively for positive and negative energy states in the I = —n — 1 
channel. The superscripts '±' means that the value of tan fi± is taken in (24). The values 
of E^jf, Q, and A are then obtained by substituting (35) in (3). It should be stressed 
out that cr a (£) can have both signs depending on the sign of the mass in a 2d Dirac 
Hamiltonian. 

6. Discussion of the results.- By extending the Krein-Friedel formula to the case of a 
singular Aharonov-Bohm (AB) potential we were able to express all physical quantities 
discussed here in terms of scattering phase shifts. We have 

• obtained closed analytical form for the density of states p a {E) of the 2d Dirac Hamil- 
tonian induced by the AB potential and predicted a resonance to occur whenever a 
bound is in the spectrum; 

• calculated spectral asymmetry cr a (£), the one-loop contribution E\^ to the effective 
energy, the induced charge Q, and the axial anomaly A (cf. [22]) of the Dirac 
operator in the AB potential for a generic self-adjoint extension; 

• pointed out that the Aharonov-Casher theorem has to be corrected for singular field 
configurations; 

• determined that the phase-shift flip only occurs for positive energies. 

In a particular case, when the bound state is absent and supersymmetry holds, our result 
for Ap a (E) are consistent with an earlier calculations of the spectral asymmetry cr a (£), 
Q, and A [13]. Relation between the axial anomaly A and the phase-shift flip has been 
discussed in [23]. Their treatment, however, differs from ours since they considered a 
complementary situation of a regular finite-flux configuration and the limit k — > 0. The 
phase-shift flip cannot take place for both positive and negative energy scattering states 
(cf. Ref. [11]). One of the arguments is that in the latter case the spectral asymmetry, 
and correspondingly induced charge and the axial anomaly all would be zero. Physically, 
a bound state can occur in the spectrum when an attractive 5-function potential is put on 
top of the AB potential [16, 24]. It this case both supersymmetry and scale invariance are 
broken, the differential scattering cross section is asymmetric, and the Hall effect occur 
[24]. 
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Note that one has the unitary equivalence between a spin 1/2 charged particle in a 
2d magnetic field and a spin 1/2 neutral particle with an anomalous magnetic moment 
in a 2d electric field [25] and our results apply to the this case as well. An application of 
our results to persistent currents in the field of a cosmic string [2, 3] and in the field of a 
gravitational vortex, and the discussion of the total energy in the AB potential are given 
elsewhere [24]. 

I should like to thank A. Comtet, Y. Georgelin, M. Knecht, S. Ouvry, and J. Stern 
for many useful and stimulating discussions, and R. C. Jones for careful reading of the 
manuscript. 
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